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ABSTRACT 

A modeling framework for the internal conformational dynamics and external 
mechanical movement of single biological macromolecules in aqueous solution 
at constant temperature is developed. Both the internal dynamics and external 
movement are stochastic; the former is represented by a master equation for a set 
of discrete states, and the latter is described by a continuous Smoluchowski equa- 
tion. Combining these two equations into one, a comprehensive theory for the 
Brownian dynamics and statistical thermodynamics of single macromolecules 
arises. This approach is shown to have wide applications. It is applied to 
protein-ligand dissociation under external force, unfolding of polyglobular pro- 
teins under extension, movement along linear tracks of motor proteins against 
load, and enzyme catalysis by single fluctuating proteins. As a generalization 
of the classic polymer theory, the dynamic equation is capable of characterizing 
a single macromolecule in aqueous solution, in probabilistic terms, (1) its ther- 
modynamic equilibrium with fluctuations, (2) transient relaxation kinetics, and 
most importantly and novel (3) nonequilibrium steady-state with heat dissipa- 
tion. A reversibility condition which guarantees an equilibrium solution and its 
thermodynamic stability is established, an iJ-theorem like inequality for irre- 
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versibility is obtained, and a rule for thermodynamic consistency in chemically 
pumped nonequilibrium steady-state is given. 

Keywords: free energy, nano-biochemistry, Smoluchowski equation, stochastic 
process, thermal fluctuation 



I. Introduction 

Progress in optics, electronics, and computer science has now made it possible to study 
biological macromolecules in aqueous solution at constant temperature by observing exper- 
imentally and measuring quantitatively the behavior of single biological macromolecules. 
These studies have been providing and will continue to yield important information on the 
behavior and properties of individual biomolecules and to reveal molecular interactions in 
and the mechanisms of biological processes. The impact which single-molecule studies will 
have on molecular biology may be gauged by comparison with the pioneering studies on 
single channel proteins in membranes, which have revolutionized physiology.^ The highly 
quantitative data obtained in these novel measurements, with piconewton and nanometer 
precision on the forces and movements of single macromolecules, complement those from 
traditional kinetic studies and structural determinations at the atomic level. 

The novel experimental approach requires a consistent theoretical framework for quanti- 
tatively understanding, interpreting, and integrating laboratory data.^^^ The objective is to 
develop a unifying molecular theory, with thermodynamic consistency, which is capable of 
integrating the three classes of quantitative measurements on macromolecules: macroscopic 
(spectroscopic kinetics), mesoscopic (single molecules), and microscopic (atomic structures). 
In this paper, we show how the spectroscopically defined kinetics, expressed in terms of dis- 
crete conformational states, is integrated with the mechanics of a macromolecule. The 
philosophy behind this approach, Stochastic Macromolecular Mechanics, is that we real- 
ize the impracticality of representing the entire conformational space of a macromolecule 
with a high-dimensional energy function. Hence we still rely on a discrete-state Markov 
model with experimentally defined "states" and kinetic parameters. However, we intro- 
duce a continuous energy landscape when there are relevant, mechanical data on the single 
molecule. Therefore the stochastic macromolecular mechanics is a mix of the discrete-state 
Markov kinetics with Brownian dynamics based on available and relevant experimental 



2 



measurements. It is a mesoscopic theory with a single set of equations. The theoretical 
approach helps researchers to identify the relevant (random) variables and key parameters 
in a macromolecular system or process, and provides them with the necessary equations. 

The discrete-state master equation approach has been long accepted as the natural 
mathematical description for biochemical kinetics of individual molecules. ^'^ With more de- 
tailed information on molecular structures and energetics, the Smoluchowski's continuous 
description of overdamped Brownian dynamics has found numerous applications in con- 
densed matter physics, polymer chemistry, and biochemistry of macromolecules^'^'^. Hanggi 
et al.^^ have reviewed the related work with Kramers' approach to chemical rate theory 
in which the assumption on overdamping is not warranted. However, for biological macro- 
molecules in aqueous solution, and with the time scale of biological interests, this assumption 
is generally acceptable. In dealing with a single protein molecule, the discrete approach is 
appropriate for spectroscopic studies^ while the continuous approach is necessary for me- 
chanical measurements. By combining these two descriptions, the stochastic macromolecu- 
lar mechanics treats the internal conformational dynamics of proteins as well as its external 
mechanics. In particular, both internal and external forces are explicitly considered. On 
the mathematical side, such a combination leads to coupled stochastic processes, ""^^ giving 
rise to three different classes of problems: reversible stationary processes (in a physicist's 
term, thermal equilibrium with fluctuations), nonstationary processes (kinetic transient), 
and irreversible stationary processes (nonequilibrium steady-state with dissipation). The 
last class of processes is novel^^"^^ and necessary for modeling motor protein (e.g. kinesin 
and myosin) movement and energetics, ^^^^^ as well as other "active macromolecules" . 

The differential equations in stochastic macromolecular mechanics are Fokker-Planck- 
master type (linear diffusion-convection equations with variable coefficients) based on con- 
servation of probability.^ This type of equations is different from the well-studied nonlinear 
reaction-diffusion equations ^^'^^ with distinctly different mathematical properties. ^^'^^ 

We would like to point out that there is already a large literature on both Smoluchowski 
equation and master equations. The novelty of our formalism is that (i) in order to com- 
bine the two types of equation into one, a condition analogous to the "thermodynamic 
box" in the elementary chemical kinetics needs to be introduced. This relation, called po- 
tential condition,^^ yields a constraint between the energy functions in the Smoluchowski 
equation and the transition rates in the master equation. This constraint guarantees the 
time-reversibility of the stationary solution to the stochastic macromolecular mechanics 
(SM3), similar to that of fluctuation-dissipation relation required in modeling equilibrium 
Brownian motion. More importantly, however, is that (ii) for single macromolecules like 
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motor proteins, this condition is violated due to the presence of chemical pumping (i.e., 
ATP hydrolysis). ^^'^^ This latter class of models based on Smoluchowski-master equation is 
relatively new. Its relation to the phenomena of stochastic resonance has been revealed only 
recently. ^^'^^ And finally, (iii) the stationary solution to the SM3 with chemical pumping 
defines a nonequilibrium steady-state the thermodynamics of which can be rigorously inves- 
tigated. By thermodynamics, we mean the entropy production, heat dissipation, nonlinear 
irreversible force-flux relationship, and the law of thermodynamics. It is this third aspect 
of the SM3, we suggest, makes the Smoluchowski approach more powerful and fundamental 
than that is generally aware. SM3 as a statistical thermodynamic theory for macromolcculcs 
in isothermal aqueous solution, both passive and active, will have wide applications, and 
deserves further investigations in the light of single-molecule experiments. 

A second objective of this paper is introducing the physical chemists who are interested 
in the Smoluchowski's approach to the unique and exciting opportunity in the current 
biophysics of nonequilibrium macromolecules on the level of single molecules. In that sense, 
SM3 is a generalization of the classic polymer theory^ into the nonlinear (Sec. II) and 
nonequilibrium (Sec. Ill) regime. 

The paper is organized with increasing complexity as follows. In Sec. II we show 
how external force is introduced into the kinetics of protein-ligand dissociation. We point 
out that mechanical measurements on single molecular complex depend critically on the 
experimental conditions - the stiffness of the force probe and the rate of its retraction. In 
Sec. Ill, polymers consisting of nonlinear subunits are introduced and we show how the 
internal kinetics is coupled to the external mechanics and movement. A Boltzmann's H- 
theorem like inequality is obtained; and the importance of nonlinear spring in serial leading 
to complex mechanical behavior is discussed. Sec. IV introduces the ATP hydrolysis into 
the model. Detailed balance, nonequilibrium steady-state, and thermodynamic consistency 
are discussed. Sec. V. shows how SM3 can be applied to the well-studied problem of 
fluctuating enzyme and yields new insights. In particular, we show how the classic concepts 
such as thermodynamic linkage and induced fit are consequence of the detailed balance, and 
can be generalized and quantified. A summary is given in Sec. VI. 



II. Macromolecular Mechanics of Protein-Ligand Dissociation 

In this section, we discuss the dissociation of a single protein-ligand complex under an 
external force introduced by an experimenter. As in any mechanics measurement, one 
first is interested in the position of the ligand with respect to the center of the mass of the 
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protein. The next mechnical quantity is what is the forces acting on the hgand. This leads 
to a Newtonian equation in which one neglects the acceleration term 

dx 

= ^intix) + Fe^X, t) +fit). (1) 

The four terms are (i) frictional force with frictional coefficient /?, {ii) intermolecular force 
between the ligand and the protein, with potential energy function Uintix): Fint{x)=—dUint{x) / dx , 
(in) the external force, and {iv) the stationary, pure random force due to collisions between 
the ligand and the solvent molecules: {f{t)) = 0. Because of the presence of the random 
force f{t), the movement x{t) is stochastic, e.g., it is a Brownian motion. Mathematically 
equivalent, the Smoluchowski's description of overdamped Brownian dynamics is based on 
a partial differential equation of parabolic type:^^'^^ 

^aPM) ^ j^^^ d^P{x,t) _ d_ j^^^^^^^^ ^ ^^^^^^^ ^^^^^^^ 

where P{x,t) is now the probability density of the ligand being at x at time t. fcg is 
the Boltzmann constant, and T is temperature which characterizes the magnitude of the 
random force f{t): {f(t)f{t')) = 2kBT6(t - t'). 

The above Eq. (|l|) and ^ lay the mathematical basis for all models, but the choices 
for Uintix) and Fg^t set the difference between different models. In the work of Shapiro 
and Qian,^^~^^ Uint = yo[{xo/xy^ — 2{xo/x)^] with a smooth repulsive force, and F^xt = 
k{x—d{t)) where k is the stiffness of the force probe exerting the external force, and d{t) = vt 
is the position of the piezoelectric motor which drives the force probe, v is the retracting 
velocity. In the work of Evans and Ritchie,^^'^'^ Unt = —Vo{xo/x)'^ {x > xq) with an abrupt 
repulsion at xq, and Fext = —F(t) is independent of x. These differences give qualitatively 
similar but qualitative different results. Hence they can be quantitatively tested against 
experimental date. 

Fig. 1 shows the results of simulations on the force-displacement curve for a protein- 
ligand complex with a simple 6-12 Lennard-Jones potential, measured using an elastic force 
probe. It is important to note that all the differences between the curves are due to the dif- 
ference in the stiffness of the force probe (k), the rate of retraction (v), and the temperature 
(ksT) at which the measurements are carried out. Therefore, this calculation demonstrates 
that the "raw" experimental data can only be understood, in general, in terms of a molec- 
ular model. It is important to realize the significance of the measurement aparatus on the 
experimental data. 



III. Macromolecular Mechanics of Polyglobular Protein Unfolding 
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In the previous section on protein- ligand dissociation, we have completely neglected the 
conformational change within the protein itself. The protein was treated as a rigid body 
exerting a force on the ligand. A more realistic model must consider the possibility of the 
protein's internal conformational change due to the external force, acting via the ligand. In 
this section, we study the unfolding of a polyglobular protein under extensional force. This 
problem naturally involves the internal dynamics of the macromolecules. 

To be concrete, let's consider the recent experimental work on giant muscle protein 
titin.'^i-^^ xitin is a protein with many globular domains (subunits) in serial. The subunits 
unfold under an external force pulling the entire molecule. The folded state of each subunit 
is rigid, and the unfolded state of each subunit can be regarded as a coiled polymer spring. 
Hence the conformational state of the entire protein, to a first order approximation, can 
be characterized by n: the number of unfolded subunits within the molecule. Let's assume 
the total number of subunits are A^, and let x be the total extension of the titin molecule 
(along the axis of external force), then a realistic characterization of a titin molecule is by 
two dynamic variables (x,n), (x > 0, < n < N). 

The equation of motion for x is again 



in which n is itself a random (discrete-state Markov) process. Hence the above equation is 
coupled to a master equation 



where P{n,t) is the probability for n at time t. Xf and are folding and unfolding rate 
constants of individual subunits. They are functions of the force acting on the subunit, 
which in turn is determined by the total extension of the molecule (x) and the number of 
unfolded domains in the chain (n). 

A comprehensive description of both the internal dynamics and external movement can 
be obtained by combining Eq. (^) and (^). We therefore have 




(3) 



dP{n,t) 
di 



{N -n + l)Xu{x, n - l)P{n - 1, t) + (n + l)A/(x, n + l)P(x, n + 1, t) 



[nXf{x, n) + (A^ — n)Xu{x, n)]P{n, t) 



(4) 



dP{x, n, t) 
di 



\ 13 J dx"^ j3 dx \ dx 

+ {N -n + l)Xu{x,n -l)P{x,n -l,t) 




+ {n + l)Xf{x,n + l)P{x,n + l,t) 



— [nAj-(x, n) + {N — n)Xu{x, n)]P{x, n, t) 



(5) 
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where P{x, n, t) is the joint probability distribution of the titin molecule having inter- 
nally n unfolded domains and external extension x. 

As in the previous section, particular models will provide specific A/(x, n), Xu{x, n), and 
Uint{x,n). These functions are not totally independent, however. Microscopic reversibility 
dictates that 



Xu{x,n) _ n+1 
Xf{x, n + 1) N — n 



exp 



Uint{x, n + 1) - Uint{x, n) 



(6) 



This condition guarantees that the stationary solution to Eq. (|5|) is a thermodynamic 
equilibrium with time reversibility. Furthermore, this reversibility condition (also known 
as potential condition and detailed balance) also guarantees the thermodynamic stability of 
the equilibrium state in terms of a generalized free energy function. As we shall see below, 
without (^) , the stationary solution in general represents a nonequilibrium steady-state with 
dissipation. 

With the reversibility condition and in the absence of external force Fext-, it is easy to 
verify that the stationary solution to Eq. (^) is 



(7) 



where 



N 
n=0 



-U{x,n)/kBT 



The time-dependent solutions to (|5|) are dynamic models for nonstationary transient kinetics 
of macromolecules. 

We now show that the equilibrium solution P*{x, n) is asymptotically stable. By stabil- 
ity, we mean a molecule approaches to its equilibrium state irrespective of its initial state. 
We introduce a free energy functional: 

^'[P(x,n,t)] = V/ {U{x,n)P{x,n,t) +kBTP{x,n,t)lnP{x,n,t))dx 



n=0' 



N .00 

-kBT\nZ + kBTy" / P{x,n,t)\n 



P{x, n, t) 
P*{x, n) 



dx 



in which the second term (known as relative entropy ^^'^^ ) is always nonnegative and equal 
to zero if and only if P{x, n, t) = P*{x, n).^^ Based on Eq. (|5|) the time derivative of ^' is 



^Pix,n,t)] = -J2 / 



I3^^j\x, t)P~^ + ksT ( J+ - J-) In 



Jn 



dx<0 (8) 
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where 

jt ,N , rT.dP{x,n,t) dU{x,n) 
J{x,t) = -ksT — — — P{x,n,t) 

and 

J+(t) = {n + l)Xf{x,n + l)P{x,n + l,t), J-{t) = {N - n)Xu{x,n)P{x,n,t). 

The integrand in Eq. (^) is always positive. Hence ^' is a Lyapunov functional for the 
time-dependent solution of Eq. (|5|), which guarantees P* to be asymptotically stable. Fur- 
thermore, because of the reversibility condition, the differential equation in (^) is symmetric 
and all its eigenvalues are real, indicating relaxations in such systems can not oscillate. 

The physical interpretation of the above result is important: It relates the equation 
of SM3 to the second law of thermodynamics. The ^'-function is the generalization of the 
equilibrium free energy of a closed, isothermal molecular system. ^ decreases monotonically 
to its minimum —ksT In Z , the Gibbs free energy, when the system reaches its equilibrium. 
It is intriguing to note that the dynamics in Eq. (^) is not governed by the gradient of the 
free energy. Nevertheless, one should see the analogy between Eq. and the H-theovem. 
of Boltzmann in his approach to irreversibility in an isolated system (microcanonical en- 
semble) . 

To analyze the stochastic dynamics of a complex macromolecule under extensional force, 
it is important first to have an essential understanding of its nonlinear mechanical property. 
A polyglobular protein model is a generalization of the classic "bead-and-spring" to non- 
linear spring.^ The protein subunits all have two energy minima while a simple Hookean 
spring has only one. This leads to fundamentally different behavior of the macromolecule. 
To illustrate this, let's apply the elementary Ohm's law for nonlinear springs in serial: the 
force on the springs are the same while the displacement is additive. For simplicity, assume 
each subunit has a potential function (and force) given in Fig. 2. This energy function has 
been used in recent work on globular protein folding kinetics. Then Fig. 3 gives a quanti- 
tative force-extension curve expected from a polyglobular protein with three subunits. As 
one can see, the most striking feature is the possibility of multiple branches of the curve 
with a given force. This is due to the combinatorics shown in Table I. 



IV. Macromolecular Mechanics of Motor Protein Movement 

With the presence of the reversibility (potential) condition, the previous model rep- 
resents a "passive" complex molecule. Without the external force, such molecules relax, 
multi-exponentially, to a thermodynamic equilibrium. They are biochemically interesting. 
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but they are not "alive" . One could argue that one of the fundamental properties of a living 
organism is the ability to convert energy among different forms (solar to electrical, electrical 
to chemical, chemical to mechanical, etc.). We now show how stochastic macromolecular 
mechanics can be used to develop models for chemomechanical energy transduction^^'^^ in 
single motor proteins.^^ In the absence of an external force, a motor protein is undergo- 
ing nonequilibrium steady-state with ATP hydrolysis and generating heat - representing a 
rudimentary form of energy metabolism. 

The key for developing a theory for motor protein is to consider that while biochemical 
studies of a protein in test tubes probe a set of discrete conformational states of the molecule, 
the mechanical studies of a protein mcasTirc positions and forces. Internally, a motor protein 
has many different conformational states within a hydrolysis cycle, and a reaction scheme 
usually can be deduced from various kinetic studies. While the protein is going through 
its conformational cycles, its center of mass moves along its designated linear track (e.g., 
kinesin on a microtubule, myosin on an actin filament, and polymerase on DNA) which 
usually has a periodic structure. The movement is stochastic; the interaction between the 
motor and the track (the force field) are usually different for different internal states of the 
molecule. 

These basic facts lead to the following equation 

dP{x,n,t) fkBT\d^P{x,n,t) ^1 d \/dUint{x,n) ^„ \ , J 

N 

+ Y.\.^kn{x)P{x,k,t)-\nk{x)P{x,n,t)\, (9) 

k=l 

where P{x, n, t) denote the joint probability of a motor protein with internal state n and 
external position x. Uintix,n) is the interaction energy between the protein in state n and 
the track. X^mix) is the transition rate constant from internal state i to state m when the 
protein is located at x. Some of the A's are pseudo-first order rate constants which contain 
the concentrations [ATP], [ADP], and [Pi]. 

When the ratio [ADP] [Pi] /[ATP] = Kgg, the equilibrium constant for the hydrolysis 
reaction, the U{x,n) and A|^(x) are again constrained by the reversibility (the superscript 
* is to indicate that the pseudo-first order rate constants are calculated in terms of the 
equilibrium concentrations) : 

Kmix) _ / Uintix,m) -Uint{xJ) \ ^^^^ 



This relation was called "thermodynamic-consistency" by T. L. Hill'^'' in his landmark con- 
tribution to the Huxley's theory of muscle contraction^^. It has to be satisfied by every 
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motor protein models. It is clear, however, that the ATP, ADP, and Pi can be kept at arbi- 
trary values. Hence in general the stationary solution of the Eq. (|9|) will be a nonequilibrium 
steady-state. Sustaining the concentrations is a form of "pumping" which keeps the system 
at nonequilibrium steady-state'^^ with positive entropy production and heat generation. ^^'^^ 
Such a molecular device is also known as isothermal ratchet. ^^'^^'^^ See Ref. 21 and 15 for 
reviews of the vast literature. If the concentrations are not actively sustained, then they 
will change slowly (since there is only a single molecule at work hydrolyzing ATP) and 
eventually reach a thermal equilibrium state in which (|l^) is satisfied. 

As in the previous sections, a practical model requires specific choices for the parameters 
U and A. In the past several years, a large amount of work have appeared on modeling 
translational motor proteins such as myosin and kinesin^^'^^'^^'^^'^^'^^ and rotational motor 
proteins such as ATP synthase. ^^'^'^ A thermodynamically valid model for a motor protein 
has to satisfy Eq. ( p!o[ ) when [ADP] [Pi] / [ATP] = K(,q, but in general without a potential 
function. This rule has not been enforced in some of the models. 

Simpler but phenomenological models based on discrete-state kinetics have also been 
developed for motor protein kinetics and energetics. ^'^'^^"■^^ It is important to point out that 
these models are completely in accord with the present continuous theory. However, drastic 
simplifications are used in order to make the models more accessible to experimental data. 
Fig. 4 shows the conceptual relationship between these two classes of models. Therefore, the 
discrete model should not be viewed as an alternative to the SM3, rather it is a simplification 
which can be further scrutinized in terms of the general theory of SM3. The stistical 
thermodynamics associated with the continuous approach, however, can be developed in 
parallel for the discrete models. 



V. Macromolecular Mechanics of Fluctuating Enzyme 

Equilibrium conformational fluctuation of proteins play an important role in enzyme 
kinetics. The theory of fluctuating enzyme^^ can be developed naturally in terms of the 
above equations for stochastic macromolecular mechanics. Let's consider a single enzyme, 
with its internal conformation characterized by x, and number of substrate molecules. 
The enzyme catalyzes a reversible isomerization reaction between two forms of the substrate 
(reactant and product), with rate constants \+{x) and A_(x). 

The equation for the catalytic reaction coupled with the enzyme conformational fluctu- 
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ation, according to stochastic macromolecular mechanics, is 
dP{n, X, t) 



dt 



-[nX+ + {N - n)A_]P(n, x, t) + (n + 1)A+P(n + 1, x, t) 



ffi Pin T f\ r) 

+ [N-n + l)A_P(n - 1, x, t) + D \ I ' ' + k-^{xP{n, x, t)), (11) 

ax^ ox 

{0<n<N) 

where P{n, x, t) is the probabihty of at time t having n number of reactant molecules and the 
enzyme internal conformation being at x. D and k characterize the protein conformational 
fluctuation, x is perpendicular to the isomerization reaction coordinate as first proposed 
by Agmon and Hopfield/'^ in contrast to the other models which address random energy 
landscape along the reaction coordinated^ Eq. (^TJ), which is essentially the same equation 
for the modeling of polyglobular protein unfolding (Eq. ||) , unifies and generalizes most of 
the important works on fluctuating enzymes. 

Along this approach, most work in the past have addressed the non-stationary, time- 
dependent solution to ([ll|). These studies are motivated by macroscopic experiments which 
are initiated (t = 0) with all the substrate in only the reactant form. If A_(x) = and 
A'^ = 1, Eq. ( |Tll) is reduced to that of Agmon and Hopfield.^'^ If A_(x) = but N is large, 
then one can introduce a continuous variable ^ = n/N, known as the survival probability, 
and Eq. (Ill]) can be approximated as (see Appendix I for more discussion) 



= A+(x)|(me,:r,t)) + D^-^§^ + k^ixPitx,t)). (12) 

Att = 0, Prob{^ = 1} = 1. 

The moments of ^, {C^){x, t) = Jq C'^P{£,, x, t)d^, can be easily obtained from Eq. (12): 

= -mX^ix){n + D-^ + k- (xin) ■ (13) 

Note Eq. (|l|) with A = 1 and Eq. (|l|) for (^) are idential. For A+(x) = ax^ + /3x + 7, 
Eq. ([1^ ) can be exactly solved by various methods if one realizes that its solution has 
a Gaussian ioim^^~^^ (also see Appendix II). From Eq. (|T3|) one immediately sees that 
high-order moments {^"^) is related to (^) by A+(x) mX^{x).^^ 

Different choices for A+(x) lead to quantitatively different models for fluctuating enzymes 
in the literature. X+{x) oc e~"^ represents a fluctuating activation energy barrier;^*' A+(x) = 
a{x+x) (> 0) representing a fluctuating cofactor concentration;^'' A+(x) = ax^ representing 
a fluctuating geometric bottleneck. 
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We now consider the reversible reaction (with A_(a;) 7^ 0) which has not been discussed 
previously. This class of models is more appropriate for recent measurements in single- 
molecule enzymology.'^ Again we assume being large. Hence we have 

dP{tx,t) d( ^ dP{^,x,t) \ , d 



dt V ^ J di 

+ Z?^!^i^ + ^(l-.(x, e)P(e, X, t)) (14) 

where D^{i,x) = + (1 - 0A-]/2A and = i\+ - (1 - 0^-- Eq. (|l|) is a 2D 
diffusion-convection equation similar to a continuous model we proposed for motor protein 
movement. One important consequence of this formulation and the reversibility condi- 
tion is realizing that conformational fluctuations of the enzyme, Vx can not be independent 
of the substrate. This constitutes the essential idea of induced fit^^~^^ and thermody- 
namic linkage!^^'^'^ For equilibrium fluctuation, again reversibility (i.e., potential condition) 
dictates:^"*^ 

dx \D^{i,x)) [^A+ + (l-0A-]2 ^ ^ 

where A' = dX{x)/dx. Therefore, 

where Vq{x) is an arbitrary function of x but it is independent of ^. As can be seen, if 
A_ << A+, then there is no requirement for ^-dependent Vx- 



VI. Conclusions 



Biological macromolecules are the cornerstone of molecular biology. Mathematical mod- 
eling of biomolecular processes requires a comprehensive and thermodynamically consistent 
theoretical basis upon which quantitative analyses can be carried out and rigorously com- 
pared with experiments. In this paper, a formal theory, we call stochastic macromolec- 
ular mechanics, is presented. The theory offers a dynamic equation for describing the 
internal kinetics as well as external motion of macromolecules in aqueous solution at con- 
stant temperature. Systematically applying this theory to various biomolecular processes 
will bring molecular biophysics closer to the standard of theoretical chemistry and physics. 



12 



At present time, Smoluchowski is well-known for its importance in calculating the micro- 
scopic fluctuations of an isothermal equilibrium system. It is less known that it can 
also be a cogent model for a macromolecules under chemical pumping. ^^'^^•'^'^'^^ What has 
not been appreciated is that this mesoscopic model also yields equilibrium and nonequi- 
librium thermodynamics for the macromolecule. Therefore, it deserves the same status as 
that of Newton's for mechanics, Navier-Stokes' for fluid dynamics. Maxwell's for electrody- 
namics, Schrodinger's for quantum mechanics, and Boltzmann's for statistical mechanics of 
isolated systems. 
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Appendix I 

Let's use the well-known linear death process^^ as an example to illustrate the continuous 
approximation for the discrete model: 

= -nAP„(t) + (n + l)AP„+i(t) (17) 

where Pn{t) is the probability of survival population being n at time t. The solution to this 
equation is well known^^ 

where is the total population at time t = 0. It is easy to show that the moments 
{nit)) = Ne-^\ 
{nHt)) 



1 + 



(n(t))2 N 



^(e--l), (18) 



(n(t))3 ^ A2 A2 \^ ■ 



We now consider the continuous counterpart of ( [T7| ) with ^ = n/N: 
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which has solution 



for initial condition 0) = — 1). The moments for ^ are 



-kXt 



(19) 



Comparing (18) and (19), we note that the continuous approximation is valid when the N 
is large and t is small. More precisely, InA^ ^ At. 



Appendix II 



Let's consider the following equation 



du 
'dt 



d 



ox^ ox 



In Gaussian form^^ which is equivalent to path integral calculation^*^ u(x, t) = ^lirait)] 
exp {^{t) — [x — ^(t)]2/2(T2(t)}, and equate coefficients of like order terms in x we have 

du/dt = -a[/i2(t) +o-2(t)] -/3|i(t) -7 
d^i/dt = -[k + 2a(j'^{t)]iJL{t) - I3(j'^{t) 
da'^/dt = 2D -2ka'^{t) -2aa^{t) 

with initial condition (t2(0) = D/k, ^(0) = 0, and z^(0) = 0. We thus have 

{u + k) + {uj- A;)e-2'^* 



2D 



(w + A;)2 - (a; - k^e 



2/3 D\ -l + e-'^^ 

) {uj + k) -{oj - k)e-'^^ 



UJ 



'P^D oj-k \ 4f3'^D 
7 M + 



u){u> — k) 



{uj + k) - {lo - k)e-'^^ 



1 

2k 



2 



4ujk 



2„-2Lut 



where u'^ = k'^ + 4Da. 
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Figure Captions 



Figure 1. The force-displacement curve for a simple Lennard- Jones bond calculated based 
on the solution of stochastic dynamics Eq. @. The Fint = — \/x^^ which equals zero 
when a; = 1, and F(.xt{x) = k{x — vt — 1). The ordinate is [F) = Fext{x)P{x,t)dx, and 
the abscissa is {x). The line labeled LJ is the expected Lennard- Jones force. As one can 
see, with increasing temperature, decreasing retracting rate, and stiffer probe, the measured 
force-displacement curve approaches to the LJ curve. It is also noted that with small ksT, 
there is a mechanical "bond rupturing" while at larger ksT, the maximal force can be 
significantly smaller than that of LJ due thermally activated transtion. 

Figure 2. (A) A conceptual energy landscape for a globular protein, cooperative folding^^"^^ 
(not to the appropriate scale) with a transition state at x^. For a single domain of titin, 
the reasonable Xo = WA and L = 30nm. The folded state is represented by the deep 
energy well and unfolded (random coil) state is represented by a shallow well with large 
entropy. The reaction coordinate (the abscissa) is uniquely defined by the direction of the 
mechanical force which pulls the molecule. At very small x, there is a closely packed core 
of all the atoms in the domain. Large x asymptotically approaches the contour length L of 
the polypeptide chain: 

ksTL 



E{x) = a ■ 



Vo 

Xq 



-2(- + 

X J \ X 



I3x + ^ 



/p(l-x/L)_ 

where Xq is the size of folded protein, Vq is the energy of the folded state, and ip is the 
persistence length of the polypeptide random coil, a, (3 and 7 are parameters characterizing 
folded, unfolded, and stretched states of the molecule. In the figure they are chosen as 
oiVq/xo = 20, (3 = 1, and jkBTL/£p = 1, and L = 10. (B) The corresponding force as 
function of displacement, -^1(2;) = dE{x) / dx. /, //, and /// are used to label the three 
monotonic regions of the curve. Note region // including the transition state is mechanically 
unstable. 

Figure 3. The force-extension relation for a series of three globular domains each of which 
is characterized by Fig. 2B. The integers by the curves are labels for the branches in the 
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curve. For three subunits in serial, the force is the same on different subunits and the total 
extension is the sum of the three individual extensions. Each branch is a sum of three from 
the /, //, and /// in Fig. 2B. Therefore, there arc total 27 branches for a trimer, but only 
10 are distinguishable (see Table I). Any branch involves // (dotted lines) is mechanically 
unstable. The stable branches are 1, 5, 8, 10. The dashed lines with the slope —k represent 
the force-displacement for the elastic force probe with stiffness k.^^'^^ The bold sawtooth 
curve^^ is expected from a mechanical force-extension experiment. A measurement using a 
force probe with less stiffness and slower rate will show less of the sawtooth pattern. 

Figure 4. x = n — 1, n, n + 1, ... in the figure represent the periodic binding sites of a 
motor protein along its track. A, B and C are the chemical states of the motor protein, 
i.e., the cyclic hydrolysis reaction can be written as A ^ B ^ C ^ A. If the potential 
energy Uint{x, k) is such that the motor can move along the track only simultaneously when 
B — > C, and there is well-defined energy barriers between x = n and x = n -|- 1 for ^, 
B, and C, then we have a simplified discrete model (line with bold face) for the stochastic 
kinetics of a motor protein. One of the most important consequences of these assumptions 
is that the ATP hydrolysis and the motor protein stepping are tightly coupled. 
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Table I 



curve # 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


composition 


1,1,1 


I,I,II 


I,II,II 


11,11,11 


I,I,III 


I,II,III 


II,II,III 


I,III,III 


II,III,III 


111,111,111 


multiplicity 


1 


3 


3 


1 


3 


6 


3 


3 


3 


1 



Table Caption 

Table I Each branch in Fig. 3 consists of a sum of 3 terms in Fig. 2B: I, //, and /// 
(Colume 2). There are total 27 branches, but some of them overlap and Colume 3 shows 
the multiplicity. 
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Figure 2: 
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